A high heterozygosity genome assembly of Aedes albopictus enables the discovery of the association of PGANT3 with blood-feeding behavior.

The Asian tiger mosquito, Aedes albopictus, is a global invasive species, notorious for its role in transmitting dangerous human arboviruses such as dengue and Chikungunya. Although hematophagous behavior is repulsive, it is an effective strategy for mosquitoes like Aedes albopictus to transmit viruses, posing a significant risk to human health. However, the fragmented nature of the Ae. albopictus genome assembly has been a significant challenge, hindering in-depth biological and genetic studies of this mosquito. In this research, we have harnessed a variety of technologies and implemented a novel strategy to create a significantly improved genome assembly for Ae. albopictus, designated as AealbF3. This assembly boasts a completeness rate of up to 98.1%, and the duplication rate has been minimized to 1.2%. Furthermore, the fragmented contigs or scaffolds of AealbF3 have been organized into three distinct chromosomes, an arrangement corroborated through syntenic plot analysis, which compared the genetic structure of Ae. albopictus with that of Ae. aegypti. Additionally, the study has revealed a phylogenetic relationship suggesting that the PGANT3 gene is implicated in the hematophagous behavior of Ae. albopictus. This involvement was preliminarily substantiated through RNA interference (RNAi) techniques and behavioral experiment. In summary, the AealbF3 genome assembly will facilitate new biological insights and intervention strategies for combating this formidable vector of disease. The innovative assembly process employed in this study could also serve as a valuable template for the assembly of genomes in other insects characterized by high levels of heterozygosity.

Targeting inflammasomes and pyroptosis in retinal diseases-molecular mechanisms and future perspectives.

Retinal diseases encompass various conditions associated with sight-threatening immune responses and are leading causes of blindness worldwide. These diseases include age-related macular degeneration, diabetic retinopathy, glaucoma and uveitis. Emerging evidence underscores the vital role of the innate immune response in retinal diseases, beyond the previously emphasized T-cell-driven processes of the adaptive immune system. In particular, pyroptosis, a newly discovered programmed cell death process involving inflammasome formation, has been implicated in the loss of membrane integrity and the release of inflammatory cytokines. Several disease-relevant animal models have provided evidence that the formation of inflammasomes and the induction of pyroptosis in innate immune cells contribute to inflammation in various retinal diseases. In this review article, we summarize current knowledge about the innate immune system and pyroptosis in retinal diseases. We also provide insights into translational targeting approaches, including novel drugs countering pyroptosis, to improve the diagnosis and treatment of retinal diseases.

De novo diploid genome assembly using long noisy reads.

The high sequencing error rate has impeded the application of long noisy reads for diploid genome assembly. Most existing assemblers failed to generate high-quality phased assemblies using long noisy reads. Here, we present PECAT, a Phased Error Correction and Assembly Tool, for reconstructing diploid genomes from long noisy reads. We design a haplotype-aware error correction method that can retain heterozygote alleles while correcting sequencing errors. We combine a corrected read SNP caller and a raw read SNP caller to further improve the identification of inconsistent overlaps in the string graph. We use a grouping method to assign reads to different haplotype groups. PECAT efficiently assembles diploid genomes using Nanopore R9, PacBio CLR or Nanopore R10 reads only. PECAT generates more contiguous haplotype-specific contigs compared to other assemblers. Especially, PECAT achieves nearly haplotype-resolved assembly on B. taurus (Bison×Simmental) using Nanopore R9 reads and phase block NG50 with 59.4/58.0 Mb for HG002 using Nanopore R10 reads.

Single-cell analysis of isoform switching and transposable element expression during preimplantation embryonic development.

Alternative splicing is an essential regulatory mechanism for development and pathogenesis. Through alternative splicing one gene can encode multiple isoforms and be translated into proteins with different functions. Therefore, this diversity is an important dimension to understand the molecular mechanism governing embryo development. Isoform expression in preimplantation embryos has been extensively investigated, leading to the discovery of new isoforms. However, the dynamics of isoform switching of different types of transcripts throughout the development remains unexplored. Here, using single-cell direct isoform sequencing in over 100 single blastomeres from the mouse oocyte to blastocyst stage, we quantified isoform expression and found that 3-prime partial transcripts lacking stop codons are highly accumulated in oocytes and zygotes. These transcripts are not transcription by-products and might play a role in maternal to zygote transition (MZT) process. Long-read sequencing also enabled us to determine the expression of transposable elements (TEs) at specific loci. In this way, we identified 3,894 TE loci that exhibited dynamic changes along the preimplantation development, likely regulating the expression of adjacent genes. Our work provides novel insights into the transcriptional regulation of early embryo development.

DNA 5-methylcytosine detection and methylation phasing using PacBio circular consensus sequencing.

Long single-molecular sequencing technologies, such as PacBio circular consensus sequencing (CCS) and nanopore sequencing, are advantageous in detecting DNA 5-methylcytosine in CpGs (5mCpGs), especially in repetitive genomic regions. However, existing methods for detecting 5mCpGs using PacBio CCS are less accurate and robust. Here, we present ccsmeth, a deep-learning method to detect DNA 5mCpGs using CCS reads. We sequence polymerase-chain-reaction treated and M.SssI-methyltransferase treated DNA of one human sample using PacBio CCS for training ccsmeth. Using long (≥10 Kb) CCS reads, ccsmeth achieves 0.90 accuracy and 0.97 Area Under the Curve on 5mCpG detection at single-molecule resolution. At the genome-wide site level, ccsmeth achieves >0.90 correlations with bisulfite sequencing and nanopore sequencing using only 10× reads. Furthermore, we develop a Nextflow pipeline, ccsmethphase, to detect haplotype-aware methylation using CCS reads, and then sequence a Chinese family trio to validate it. ccsmeth and ccsmethphase can be robust and accurate tools for detecting DNA 5-methylcytosines.

High-throughput and high-accuracy single-cell RNA isoform analysis using PacBio circular consensus sequencing.

Although long-read single-cell RNA isoform sequencing (scISO-Seq) can reveal alternative RNA splicing in individual cells, it suffers from a low read throughput. Here, we introduce HIT-scISOseq, a method that removes most artifact cDNAs and concatenates multiple cDNAs for PacBio circular consensus sequencing (CCS) to achieve high-throughput and high-accuracy single-cell RNA isoform sequencing. HIT-scISOseq can yield >10 million high-accuracy long-reads in a single PacBio Sequel II SMRT Cell 8M. We also report the development of scISA-Tools that demultiplex HIT-scISOseq concatenated reads into single-cell cDNA reads with >99.99% accuracy and specificity. We apply HIT-scISOseq to characterize the transcriptomes of 3375 corneal limbus cells and reveal cell-type-specific isoform expression in them. HIT-scISOseq is a high-throughput, high-accuracy, technically accessible method and it can accelerate the burgeoning field of long-read single-cell transcriptomics.

High-throughput Pore-C reveals the single-allele topology and cell type-specificity of 3D genome folding.

Canonical three-dimensional (3D) genome structures represent the ensemble average of pairwise chromatin interactions but not the single-allele topologies in populations of cells. Recently developed Pore-C can capture multiway chromatin contacts that reflect regional topologies of single chromosomes. By carrying out high-throughput Pore-C, we reveal extensive but regionally restricted clusters of single-allele topologies that aggregate into canonical 3D genome structures in two human cell types. We show that fragments in multi-contact reads generally coexist in the same TAD. In contrast, a concurrent significant proportion of multi-contact reads span multiple compartments of the same chromatin type over megabase distances. Synergistic chromatin looping between multiple sites in multi-contact reads is rare compared to pairwise interactions. Interestingly, the single-allele topology clusters are cell type-specific even inside highly conserved TADs in different types of cells. In summary, HiPore-C enables global characterization of single-allele topologies at an unprecedented depth to reveal elusive genome folding principles.

NanoSNP: a progressive and haplotype-aware SNP caller on low-coverage nanopore sequencing data.

MOTIVATION: Oxford Nanopore sequencing has great potential and advantages in population-scale studies. Due to the cost of sequencing, the depth of whole-genome sequencing for per individual sample must be small. However, the existing single nucleotide polymorphism (SNP) callers are aimed at high-coverage Nanopore sequencing reads. Detecting the SNP variants on low-coverage Nanopore sequencing data is still a challenging problem. RESULTS: We developed a novel deep learning-based SNP calling method, NanoSNP, to identify the SNP sites (excluding short indels) based on low-coverage Nanopore sequencing reads. In this method, we design a multi-step, multi-scale and haplotype-aware SNP detection pipeline. First, the pileup model in NanoSNP utilizes the naive pileup feature to predict a subset of SNP sites with a Bi-long short-term memory (LSTM) network. These SNP sites are phased and used to divide the low-coverage Nanopore reads into different haplotypes. Finally, the long-range haplotype feature and short-range pileup feature are extracted from each haplotype. The haplotype model combines two features and predicts the genotype for the candidate site using a Bi-LSTM network. To evaluate the performance of NanoSNP, we compared NanoSNP with Clair, Clair3, Pepper-DeepVariant and NanoCaller on the low-coverage (∼16×) Nanopore sequencing reads. We also performed cross-genome testing on six human genomes HG002-HG007, respectively. Comprehensive experiments demonstrate that NanoSNP outperforms Clair, Pepper-DeepVariant and NanoCaller in identifying SNPs on low-coverage Nanopore sequencing data, including the difficult-to-map regions and major histocompatibility complex regions in the human genome. NanoSNP is comparable to Clair3 when the coverage exceeds 16×. AVAILABILITY AND IMPLEMENTATION: https://github.com/huangnengCSU/NanoSNP.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Applications and potentials of nanopore sequencing in the (epi)genome and (epi)transcriptome era.

The Human Genome Project opened an era of (epi)genomic research, and also provided a platform for the development of new sequencing technologies. During and after the project, several sequencing technologies continue to dominate nucleic acid sequencing markets. Currently, Illumina (short-read), PacBio (long-read), and Oxford Nanopore (long-read) are the most popular sequencing technologies. Unlike PacBio or the popular short-read sequencers before it, which, as examples of the second or so-called Next-Generation Sequencing platforms, need to synthesize when sequencing, nanopore technology directly sequences native DNA and RNA molecules. Nanopore sequencing, therefore, avoids converting mRNA into cDNA molecules, which not only allows for the sequencing of extremely long native DNA and full-length RNA molecules but also document modifications that have been made to those native DNA or RNA bases. In this review on direct DNA sequencing and direct RNA sequencing using Oxford Nanopore technology, we focus on their development and application achievements, discussing their challenges and future perspective. We also address the problems researchers may encounter applying these approaches in their research topics, and how to resolve them.

Environmental Water and Sediment Microbial Communities Shape Intestine Microbiota for Host Health: The Central Dogma in an Anthropogenic Aquaculture Ecosystem.

From increasing evidence has emerged a tight link among the environment, intestine microbiota, and host health status; moreover, the microbial interaction in different habitats is crucial for ecosystems. However, how the environmental microbial community assembly governs the intestinal microbiota and microbial communities of multiple habitats contribute to the metacommunity remain elusive. Here, we designed two delicate experiments from temporal and spatial scales in a shrimp culture pond ecosystem (SCPE). Of the SCPE metacommunity, the microbial diversity was mainly contributed to by the diversity of-ß IntraHabitats and ß InterHabitats , and water and sediment communities had a large contribution to the shrimp intestine community as shown by SourceTracker and Sloan neutral community model analyses. Also, phylogenetic bin-based null model results show that microbial assembly of three habitats in the SCPE appeared to be largely driven by stochastic processes. These results enrich our understanding of the environment-intestinal microbiota-host health closely linked relationship, making it possible to be the central dogma for an anthropogenic aquaculture ecosystem. Our findings enhance the mechanistic understanding of microbial assembly in the SCPE for further analyzing metacommunities, which has important implications for microbial ecology and animal health.

Genome-wide detection of cytosine methylations in plant from Nanopore data using deep learning.

In plants, cytosine DNA methylations (5mCs) can happen in three sequence contexts as CpG, CHG, and CHH (where H = A, C, or T), which play different roles in the regulation of biological processes. Although long Nanopore reads are advantageous in the detection of 5mCs comparing to short-read bisulfite sequencing, existing methods can only detect 5mCs in the CpG context, which limits their application in plants. Here, we develop DeepSignal-plant, a deep learning tool to detect genome-wide 5mCs of all three contexts in plants from Nanopore reads. We sequence Arabidopsis thaliana and Oryza sativa using both Nanopore and bisulfite sequencing. We develop a denoising process for training models, which enables DeepSignal-plant to achieve high correlations with bisulfite sequencing for 5mC detection in all three contexts. Furthermore, DeepSignal-plant can profile more 5mC sites, which will help to provide a more complete understanding of epigenetic mechanisms of different biological processes.

SCSit: A high-efficiency preprocessing tool for single-cell sequencing data from SPLiT-seq.

SPLiT-seq provides a low-cost platform to generate single-cell data by labeling the cellular origin of RNA through four rounds of combinatorial barcoding. However, an automatic and rapid method for preprocessing and classifying single-cell sequencing (SCS) data from SPLiT-seq, which directly identified and labeled combinatorial barcoding reads and distinguished special cell sequencing data, is currently lacking. Here, we develop a high-efficiency preprocessing tool for single-cell sequencing data from SPLiT-seq (SCSit), which can directly identify combinatorial barcodes and UMI of cell types and obtain more labeled reads, and remarkably enhance the retained data from SCS due to the exact alignment of insertion and deletion. Compared with the original method used in SPLiT-seq, the consistency of identified reads from SCSit increases to 97%, and mapped reads are twice than the original. Furthermore, the runtime of SCSit is less than 10% of the original. It can accurately and rapidly analyze SPLiT-seq raw data and obtain labeled reads, as well as effectively improve the single-cell data from SPLiT-seq platform. The data and source of SCSit are available on the GitHub website https://github.com/shang-qian/SCSit.

Integrated single-cell analysis revealed immune dynamics during Ad5-nCoV immunization.

Coronavirus disease 2019 (COVID-19), driven by SARS-CoV-2, is a severe infectious disease that has become a global health threat. Vaccines are among the most effective public health tools for combating COVID-19. Immune status is critical for evaluating the safety and response to the vaccine, however, the evolution of the immune response during immunization remains poorly understood. Single-cell RNA sequencing (scRNA-seq) represents a powerful tool for dissecting multicellular behavior and discovering therapeutic antibodies. Herein, by performing scRNA/V(D)J-seq on peripheral blood mononuclear cells from four COVID-19 vaccine trial participants longitudinally during immunization, we revealed enhanced cellular immunity with concerted and cell type-specific IFN responses as well as boosted humoral immunity with SARS-CoV-2-specific antibodies. Based on the CDR3 sequence and germline enrichment, we were able to identify several potential binding antibodies. We synthesized, expressed and tested 21 clones from the identified lineages. Among them, one monoclonal antibody (P3V6-1) exhibited relatively high affinity with the extracellular domain of Spike protein, which might be a promising therapeutic reagent for COVID-19. Overall, our findings provide insights for assessing vaccine through the novel scRNA/V(D)J-seq approach, which might facilitate the development of more potent, durable and safe prophylactic vaccines.

Genomic Elucidation of a COVID-19 Resurgence and Local Transmission of SARS-CoV-2 in Guangzhou, China.

While China experienced a peak and decline in coronavirus disease 2019 (COVID-19) cases at the start of 2020, regional outbreaks continuously emerged in subsequent months. Resurgences of COVID-19 have also been observed in many other countries. In Guangzhou, China, a small outbreak, involving less than 100 residents, emerged in March and April 2020, and comprehensive and near-real-time genomic surveillance of SARS-CoV-2 was conducted. When the numbers of confirmed cases among overseas travelers increased, public health measures were enhanced by shifting from self-quarantine to central quarantine and SARS-CoV-2 testing for all overseas travelers. In an analysis of 109 imported cases, we found diverse viral variants distributed in the global viral phylogeny, which were frequently shared within households but not among passengers on the same flight. In contrast to the viral diversity of imported cases, local transmission was predominately attributed to two specific variants imported from Africa, including local cases that reported no direct or indirect contact with imported cases. The introduction events of the virus were identified or deduced before the enhanced measures were taken. These results show the interventions were effective in containing the spread of SARS-CoV-2, and they rule out the possibility of cryptic transmission of viral variants from the first wave in January and February 2020. Our study provides evidence and emphasizes the importance of controls for overseas travelers in the context of the pandemic and exemplifies how viral genomic data can facilitate COVID-19 surveillance and inform public health mitigation strategies.

Compacting a synthetic yeast chromosome arm.

BACKGROUND: Redundancy is a common feature of genomes, presumably to ensure robust growth under different and changing conditions. Genome compaction, removing sequences nonessential for given conditions, provides a novel way to understand the core principles of life. The synthetic chromosome rearrangement and modification by loxP-mediated evolution (SCRaMbLE) system is a unique feature implanted in the synthetic yeast genome (Sc2.0), which is proposed as an effective tool for genome minimization. As the Sc2.0 project is nearing its completion, we have begun to explore the application of the SCRaMbLE system in genome compaction. RESULTS: We develop a method termed SCRaMbLE-based genome compaction (SGC) and demonstrate that a synthetic chromosome arm (synXIIL) can be efficiently reduced. The pre-introduced episomal essential gene array significantly enhances the compacting ability of SGC, not only by enabling the deletion of nonessential genes located in essential gene containing loxPsym units but also by allowing more chromosomal sequences to be removed in a single SGC process. Further compaction is achieved through iterative SGC, revealing that at least 39 out of 65 nonessential genes in synXIIL can be removed collectively without affecting cell viability at 30 °C in rich medium. Approximately 40% of the synthetic sequence, encoding 28 genes, is found to be dispensable for cell growth at 30 °C in rich medium and several genes whose functions are needed under specified conditions are identified. CONCLUSIONS: We develop iterative SGC with the aid of eArray as a generic yet effective tool to compact the synthetic yeast genome.

Efficient assembly of nanopore reads via highly accurate and intact error correction.

Long nanopore reads are advantageous in de novo genome assembly. However, nanopore reads usually have broad error distribution and high-error-rate subsequences. Existing error correction tools cannot correct nanopore reads efficiently and effectively. Most methods trim high-error-rate subsequences during error correction, which reduces both the length of the reads and contiguity of the final assembly. Here, we develop an error correction, and de novo assembly tool designed to overcome complex errors in nanopore reads. We propose an adaptive read selection and two-step progressive method to quickly correct nanopore reads to high accuracy. We introduce a two-stage assembler to utilize the full length of nanopore reads. Our tool achieves superior performance in both error correction and de novo assembling nanopore reads. It requires only 8122 hours to assemble a 35X coverage human genome and achieves a 2.47-fold improvement in NG50. Furthermore, our assembly of the human WERI cell line shows an NG50 of 22 Mbp. The high-quality assembly of nanopore reads can significantly reduce false positives in structure variation detection.

A human circulating immune cell landscape in aging and COVID-19.

Age-associated changes in immune cells have been linked to an increased risk for infection. However, a global and detailed characterization of the changes that human circulating immune cells undergo with age is lacking. Here, we combined scRNA-seq, mass cytometry and scATAC-seq to compare immune cell types in peripheral blood collected from young and old subjects and patients with COVID-19. We found that the immune cell landscape was reprogrammed with age and was characterized by T cell polarization from naive and memory cells to effector, cytotoxic, exhausted and regulatory cells, along with increased late natural killer cells, age-associated B cells, inflammatory monocytes and age-associated dendritic cells. In addition, the expression of genes, which were implicated in coronavirus susceptibility, was upregulated in a cell subtype-specific manner with age. Notably, COVID-19 promoted age-induced immune cell polarization and gene expression related to inflammation and cellular senescence. Therefore, these findings suggest that a dysregulated immune system and increased gene expression associated with SARS-CoV-2 susceptibility may at least partially account for COVID-19 vulnerability in the elderly.

Erratum: Author Correction: Immune cell profiling of COVID-19 patients in the recovery stage by single-cell sequencing.

[This corrects the article DOI: 10.1038/s41421-020-0168-9.].

Immune cell profiling of COVID-19 patients in the recovery stage by single-cell sequencing.

COVID-19, caused by SARS-CoV-2, has recently affected over 1,200,000 people and killed more than 60,000. The key immune cell subsets change and their states during the course of COVID-19 remain unclear. We sought to comprehensively characterize the transcriptional changes in peripheral blood mononuclear cells during the recovery stage of COVID-19 by single-cell RNA sequencing technique. It was found that T cells decreased remarkably, whereas monocytes increased in patients in the early recovery stage (ERS) of COVID-19. There was an increased ratio of classical CD14++ monocytes with high inflammatory gene expression as well as a greater abundance of CD14++IL1ß+ monocytes in the ERS. CD4+ T cells and CD8+ T cells decreased significantly and expressed high levels of inflammatory genes in the ERS. Among the B cells, the plasma cells increased remarkably, whereas the naïve B cells decreased. Several novel B cell-receptor (BCR) changes were identified, such as IGHV3-23 and IGHV3-7, and isotypes (IGHV3-15, IGHV3-30, and IGKV3-11) previously used for virus vaccine development were confirmed. The strongest pairing frequencies, IGHV3-23-IGHJ4, indicated a monoclonal state associated with SARS-CoV-2 specificity, which had not been reported yet. Furthermore, integrated analysis predicted that IL-1ß and M-CSF may be novel candidate target genes for inflammatory storm and that TNFSF13, IL-18, IL-2, and IL-4 may be beneficial for the recovery of COVID-19 patients. Our study provides the first evidence of an inflammatory immune signature in the ERS, suggesting COVID-19 patients are still vulnerable after hospital discharge. Identification of novel BCR signaling may lead to the development of vaccines and antibodies for the treatment of COVID-19.